UCONN-97-02 



Gauss's law, gauge invariance, and long-range forces in QCD 



Mario Belloni*, Lusheng Chen^, and Kurt Haller* 

Department of Physics, University of Connecticut, Storrs, Connecticut 06269 

(February 1, 2008) 



Abstract 



We use a unitary operator constructed in earlier work to transform the Hamil- 
tonian for QCD in the temporal {Aq = 0) gauge into a representation in which 
the quark field is gauge- invariant, and its elementary excitations — quark and 
antiquark creation and annihilation operators — implement Gauss's law. In 
that representation, the interactions between gauge-dependent parts of the 
gauge field and the spinor (quark) field have been transformed away and re- 
placed by long-range non-local interactions of quark color charge densities. 
These long-range interactions connect SU(3) color charge densities through 
an infinite chain of gauge-invariant gauge fields either to other SU(3) color 
charge densities, or to a gluon "anchor". We discuss possible implications of 
this formalism for low-energy processes, including confinement of quarks that 
are not in color singlet configurations. 
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In earlier work |],[|, we showed how to construct states that implement the non-Abelian 
Gauss's law in QCD and Yang-Mills theory, and we constructed gauge-invariant spinor 
(quark) and gauge (gluon) operator-valued fields. In the work presented here, we explore 
the implications of these formal developments for QCD dynamics. In Ref. [[| we discussed 
two unitarily equivalent representations of QCD. In one — the C representation — the 
Gauss's law operator in the temporal (Ao = 0) gauge has the form 

G a (r) = d t Ut(r)+gr bc Al(r)Wr(r)+f (r) , (1) 

where II"(r) is the momentum canonically conjugate to Af(r) as well as the chromoelectric 
field; Jq (r) is the quark color charge density j$ (r) = g ^(r) (A a /2) ?/>(r), and the A a represent 
the Gell-Mann SU(3) matrices. In the unitarily equivalent M representation, that same 
Gauss's law operator takes the form 

g a (r) = A a n?(r) = $n»(r) + gf abc A\(v) IIJ(r) . (2) 



The unitary equivalence that relates Q a (r) and Q a (r), and that accounts for the fact 
that j'o(r) does not appear explicitly in the non-Abelian Gauss's law operator in the M 
representation, is similar to a relationship between two forms of the Gauss's law operator 
in QED [|||. In QED, 0(r) = d^v) + j (r), where j (r) = e^(r)ip(r), and <?(r) is uni- 
tarily equivalent to diUi(r), which represents the complete Gauss's law operator in what 
corresponds to the Af representation for QED. When all operators and states for QED in 
the temporal gauge are subjected to the unitary transformation that transforms Q(r) into 
<9jITj(r), the unitary operator U = exp[ — i J <9jIL(r) uj(r) ] dr implements the gauge transfor- 
mations U A{(r) U^ 1 = Ai(r) + di0o(r), where oo(r) commutes with all operator-valued fields; 
and, because djli(r) is the complete Gauss's law operator in the M representation, and ^(r) 
commutes with <9jilj(r), ip(r) represents a gauge-invariant charged particle (electron) field 
in the M representation. 

When the Hamiltonian for QED in the temporal gauge is transformed into the M rep- 
resentation, the only remaining dynamical interactions between the electron field and the 
electromagnetic field are the ones between the electrons and the gauge-invariant excitations 
of the gauge field — which, in QED, correspond to transversely polarized photons 0,f§. 
The interactions which, in the original C representation, were mediated by the longitudinal, 
gauge-dependent parts of the gauge field — the interactions mediated by the exchange of 
longitudinal photon "ghosts" — appear in the M representation as the non-local Coulomb 
interaction given by 

H c = l - f J ^P^-drdr'. (3) 
2 J 4vr r-r' V ; 



The purpose of this present work is to extend this approach to the non-Abelian QCD. 
There are, of course, important differences between QED and QCD. As was shown in Refs. 
@ and 0, the states that implement Gauss's Law in the M representation for QED, i.e. 
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the states that satisfy 3jIIj(r)|n) = - or, more properly for QED, diH\ + \r)\n} = 0, 
where dill\ + ^ selects the annihilation operator part of c^IIj — constitute a Fock space that 
can easily be constructed. The A/"-representation states that satisfy Q a {r)\u) = in QCD, 
however, can only be represented as complicated coherent superpositions of perturbative 
states that are neither readily orthogonalizable nor normalizable Qj. Nevertheless, there are 
also significant similarities between QCD and QED. In particular, since Q a (r) and Q a {r) 
are unitarily equivalent, and since the quark field commutes with Af(r) and II" (r) and 
therefore with G a (r), ip(r) represents the gauge- invariant quark field in the M representation 
of QCD. And, it is possible to use the unitary equivalence demonstrated in Ref. ||, 

i a (v)=U c 1 C(r)Uc, (4) 

where £ a represents some functional of operator- valued quark and/or gluon fields in the 
C representation and £ a its transform in the M representation, to transform dynamical 
variables from one representation to the other. In Ref. |J , we showed that Q a {r) is the M 
representation form of the Gauss's law operator which, in the C representation, is given by 
Q a {r) . Here, we will apply this transformation to transform the QCD Hamiltonian from 
the C to the M representation. In the M representation, the excitations of the "bare" quark 
field ipiy) are gauge- invariant, and therefore already include the gauge-dependent gauge field 
components required to obey Gauss's law. The interactions between these components of 
the gauge fields and quarks therefore will no longer appear explicitly in the Hamiltonian in 
the M representation, but will be replaced by non-local interactions among quarks, which 
form a kind of non-Abelian analog to the Coulomb interaction in QED. 

We observe that Uq is given by Uq = e c °e c , where Co and C are 

C = i fdrX a (r) j*(r) and C = % f dry^(r) j£(r) (5) 

respectively, where X a (r) = [ j^Af(r) ] and y a (r) = ^sA?(t), and where Af(r) is defined in 
Ref. p| by an operator differential equation. We will not repeat that definition here since, 
for our purposes in this work, the form of Af(r) is less important than the fact that the 
unitary equivalence expressed in Eq. @ obtains. 

The form of the temporal gauge QCD Hamiltonian in the C representation is 

H = Jdr{ ±n?(r)n»(r) + |i^(r)^(r) 

+ ^(r) [pm-i ai (di-igA?(T)% )]^(r)} (6) 

where i^-(r) = djA$(r) - diA](r) - gf^A^Afa), and where the term J Ag(r) G a {r) dr, 
which appears in the Hamiltonian when primary constraints are avoided in the canonical 
quantization of this model ||, has been omitted because it is irrelevant for this work. It 
is obvious that, in transforming H to the M representation, F^(r) remains untransformed. 
It is much more difficult, however, to transform II" (r) to its Af representation transform, 
n"(r). We have found that IT" (r) is given by 
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n»(r) = n?(r) + E EE ^H)^/ftiSo#<)WH)«l C#)OOtfOO 



m=0 n=0 r=0 



(7) 



where f?^ 1 is the chain of SU(3) structure functions 



(v) 

j&fh _ p[\]l3b[\] y6[l]a[2]6[2] y6[2]a[3]6[3] . . . fb[ v -2]a[ v -l]b{ v -l] j%-l]a[»7]7 ? (g) 



^(r) = IlJU * aH (r) , and A^r) = Ul=i ^(r) . ^ ) (r) Jo a (r) is defined by 

^)(r)jo(r) = > (^? (a) (r) ^ (• • • «W r ) (#('))))) , (9) 

where -(r) represents the gauge- invariant gauge field, which can be expressed as 

[4u(r) $ ] = Vfc(r) [ A 6 W $ ] V c \r) + J lfc(r) ^(r) ; (10) 

Vc(r) is a unitary operator given by Vc(r) = exp ( — ^^"(rjy ) exp f— igAf'^r)^). We 
have been able to establish Eq. (0) through an inductive proof that proceeds on very similar 
lines as the proof of the "fundamental theorem" in Ref . . We have further verified Eq. (^) 
by showing that it is consistent with 

fitfl?(r) + gf ahc A 6 (r) fl?(r) + j a (r) = fiW?(r) + <7/ a6c 4(r) IT?(r) . (11) 
Equation (|7|) is a non-Abelian analog to a relation that, in QED, has the form HIJ 

n i (r) = n j (r) + ^ Jo (r), (12) 

In Eq. (|T2"D, jo( r ) is the entire electric charge density — a consequence of the fact that, in 
QED, charge resides only in the fields coupled to the gauge field, and never in the gauge 
field itself. In Eq. ([7J), jo(r) represents contributions to the SU(3) color charge density from 
quarks only. The gluon color charge density, Jq(t), does not appear explicitly in IT"(r) (or, 
later in this paper, in H), because the unitary transformation we have used to relate II" (r) 
and II"(r) does not eliminate the coupling of the gauge- dependent parts of the gauge field 
toJ a (r). 

When we transform the quark field contribution to the Hamiltonian in the M represen- 
tation, we obtain 

H quark = ft (r) (3m - iaci ( d t - ig Af(r) ^ ) ^(r) = ^ f (r) [(3m - iaA } ip(r) 

-ift(v)a t \Vc{r)d t Vc\r)-igV c (r) (4(r)|)^W]#) ; (13) 



and, from 

^(r) = V c \r) ip(r) and ^(r) = ft(r) V c (r) , (14) 

and from Eqs. (|i0|) and fll3l) , 
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H qU ark = ^ f (r) [(3m - ion (<% - W^ai( r )^)} VK r ) ■ 



(15) 



After the transformations from the C representation to the JV representation have been 
carried out on all constituent parts of H, we obtain 



H = J dr { |n?(r)IIJ(r) + \F${t)F%{t) + ^(r) [(3m - ia, (d { - igA a GH (r)%)} ^(r) 

+ 1 E EE ^(-l) m+n+ 7S/(^/ffn»(r)^f m) (r) A^(r)f| (^(r^r 



m=0 n=0 r=0 



+IEEE 

m=0 n=0 r=0 



+IEE / +r '(-l) r+r '4f/J'fl Wototfto) | (^)(r)i fe '(r))} - (16) 



r=0 r'=0 

In deriving Eq. (|IB), we have used the identity 



Eg m /_i\m npac^pp, , \ g m ' ( _i\m' rjl'ac' ~pfl' / 



m=0 



m'=0 



E E ml 



M=0m=0 



(AT-m)! 



~ L ) J(M) /C (M)' 



Since 0(M) 
reduces to 



Af M! 



m=0 m!(M-m)! 



(-!)"» = (1 - 1) M = for M > 0, while 0(0) 



E£(-i)X^W E Sr(-i) m /&o ^W r ) = * 



m=0 



'(m) ^(mjV 1 -; 

m'=0 



(17) 
1 , Eq. (0) 

(18) 



Equation (p~8|) also obtains when ■M^(r) is substituted for 72,^ (r). 

We can use the Baker-Hausdorff-Campbell (BHC) theorem to show that 



V c \r)^V c (r) = J2J2^(-l) m+ V^f^nf m) (r)Mf n) (r) 
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(19) 



m=0 n=0 



and, from Eqs. fllCf ) and ([Tj|) 



E E ^(-i) m+n /5/??^ «, jwf B) ) v? 



m=0 n=0 



m=0 n=0 



rt at rarb g m + n / -,\m+nffibcrvcd<T?fi \av \rd 

~9 A i J 1^ 2^,l*mT\~ l ) J(m) /(») K {m) ^(n) V i ■ 
m=0 n=0 

And we can use Eqs. ( |T9"D and to reexpress Eq. ( |T6| ) in the form 
H = Jdv{ |n«(r) H?(r) + |^(r)^(r) + ^(r) _ ia . - ^ H (r)f 



(20) 



V(r)} 



(21) 
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with 



dr < — Tr 



E/(- 1 ) r ^ a ( r )T^- 1 W^(r)ft(7^(r)j*(r; 



Lr=0 



-Tr 



r=0 



1^ 

2 



(22) 



and 



H 



LR 



dr { +Tr 



E? r+1 (-l)7ffrt(r)f ^- l (r)^V c (rK 14 (r)|t (7&(r)#(r) 



r=0 



+Tr 



E/ +1 (-l)7SV dCTe § (^)(r)j "(r)) K c - X (r)¥ V c (r)^ u (r)IlJ(r)^ 



Lr=0 



+ I E E / +r '(-l) r+r 7§V?f | (^,(r)j fc (r)) ft (T&Wjftr)) . (23) 



r=0 r'=0 



It is instructive to examine the physical significance of the individual terms in Eq. (pT|). 
We observe that H includes the kinetic energy of the quark and gluon fields, as well as 
an interaction term between the spatial SU(3) current term for the spinor field (which 
is gauge-invariant in the M representation) and the gauge- invariant gauge field Aq H (t), 
which is gauge- invariant in both the C and the M representations. H also includes Hg, in 
which the Gauss's law operator Q appears on either the extreme left-hand or right-hand 
side of an operator product. Hg vanishes when evaluated between states that implement 
Gauss's law in the M representation. H also includes a long-range non-local interaction, 
Hm, which consists of three terms; the last of these three terms in Eq. (p3|) describes the 
interaction between two SU(3) color charge densities through a non-local interaction that 
can be characterized as a non-Abelian analog of the Coulomb interaction in QED; but it 
differs from the Coulomb interaction in QED in that the factor (1/2) d~ 2 does not connect 
two j'o terms directly. The long-range interaction in this case is transmitted between the 
two quark color charge densities through an infinite chain of gauge-invariant gauge fields 
as shown in Eq. (^). The other two terms in Eq. (|23| ) each connect a jfi term through an 
identical chain of gauge-invariant gauge fields to the gluon field itself, "anchoring" the j'q 
term, finally, in the gluon field through an Agj(r) n|(r) term. 

It is appealing to speculate that the three terms in Eq. (|23D have an important role 
in confining quark states that are not in singlet color configurations. To be sure, such an 
idea is speculative. The spatial dependence of the long-range interactions in H LR cannot be 
determined without more information about matrix elements of gauge fields than we have 
presented here. Partial insight into these long-range interactions can be obtained by noting 
that in the virtually identical SU(2) case, the summation over all orders, s, of T£\(r) can be 
performed exactly to give 

E 9 r (-l) r 4r) = - V - Th e5dH T& + j^e s ^ db e bS ^ h T 5(1) T <5(2) , (24) 

r=0 
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where T 5 (r) = ^.(r) % and W = g 2 (T^) . W is a differential operator in which the 
dj/d 2 in each T 13 operates on all operators to its right in the series expansion of 1/(1 + W), 
as shown in Eq. (H). Equation ([24]) therefore describes a differential operator rather than 
an algebraic quantity. Nevertheless, it is a reasonable conjecture that might develop 

expectation values in the gluon "sea" in which quarks are embedded that have the effect of 
turning 1/(1 + W) into a confining potential, and that a very similar phenomenon might 
occur in the SU(3) case. We consider this scenario to be an interesting conjecture that 
warrants further study and development. 



It is of particular interest to consider Eq. ( p23|) in conjunction with a model suitable for 
low-energy QCD, that combines gluons and "static" quarks — i.e. a model in which quark 
kinetic energy is small, and quark-antiquark creation and annihilation is quenched — but 
in which Gauss's law is obeyed. This model is of interest because it allows us to exploit a 
pertinent analogy between QCD and QED. In QED, the covariant gauge formulation enjoys 
a clear advantage for perturbative calculations — particularly when radiative corrections 
are required — but is not as well adapted for extracting the Coulomb interaction as the 
dominant one for low-energy phenomena, such as the energy levels of the bound states of 
the Hydrogen atom. In contrast, in a formulation of QED in which Gauss's law has been 
implemented and the "pure gauge" components of the gauge field have been transformed 
away, the dominance of the Coulomb interaction for low-energy phenomena becomes very 
explicit Similarly, in QCD, when we eliminate the interaction between the gauge- 

dependent parts of the gauge fields and the quarks, the long-range interactions of quark 
color charge densities described in Eq. (^) become very explicit. As in the case of QED, 
these long-range interactions are not manifest in either the C-representation Hamiltonian 
described in Eq. (^|), or in the covariant-gauge formulation of QCD. It is very natural to 
suspect that these long-range interactions are important for low-energy processes, and that 
they have a connection with the color confinement observed in low-energy phenomena. A 
number of authors have suggested that there is some connection between gauge invariance 
and color confinement in non-Abelian gauge theories |B|-r2|. The following observations 



about the application of Hlr, the long-range part of the transformed Hamiltonian H, are 
pertinent to this question. 

We consider the expectation value (A| j'q (r) |A) , where |A) is an n-particle state of 
"static" quarks, in the sense defined above. We also assume that the quarks and anti- 
quarks are described by a set of wave functions w"(r) and v "(r) respectively, and that the 
quark field ^(r) is represented as 

V>(r)=£ [«(r) + gtX(r)] , (25) 

k 

where q" and q]^ annihilate quarks and create antiquarks respectively, and w"(r) and f"(r) 
are the positive and negative energy solutions of a Dirac equation in an external potential (a 
labels the SU(3) index and n the energy, angular momentum, etc. of quarks and antiquarks 
in the corresponding orbital states). We will identify |A) as a multiquark state \q"- ■ -q^,) . 
We observe that for such a state, (q%- ■ 1 / dr j'o(r) |<z"- • = when \q%- ■ -q°,) is a color 
singlet state, since / drj^r) is the generator of infinitesimal rotations in the SU(3) space, 
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and the singlet state is the scalar in that space. When 2gj(r) imposes only minimal variations 
on Jq (r) over the region in r in which u°(r) and f"(r) are localized, then the expectation 
values that enter into the evaluation of long-range forces between "static" quark states in 
color-singlet configurations can be expected to vanish. For example, when two singlet states, 
| A) = \q"- ■ -q®,) and |B) = \Q%- ■ -Qn 1 ) interact through the last term in Eq. Q23|), where the 
q% annihilate quarks centered about point R, and the annihilate quarks centered about 
R', and where |R — R'| ^> \5r\ where |d>r| is the mean radius of the wave functions w"(r) 
and v%(r), then the following conjecture becomes plausible: Because the distance between 
the end points of the chain of gauge-invariant gauge fields in Eq. (p3|) is much larger than 
l^ r l ) /(s)' l ^(s) ( r )io ( r ) wou ld terminate essentially in k j drj^r) , where k is some constant. 
In that case, the matrix element f^ h T^(r)(q"- ■ -q%, \ jft (r) \q%- ■ -q%,) and the corresponding 
chain terminating in the quark state \Q%- • w Q%) both would vanish. When |R— R'| decreases, 
the cancellation can become less complete, and "van der Waals" forces between singlet states 
become a reasonable conjecture. 

This account of how confining long-range forces might develop between multiquark states 
that are not in color singlet configurations, is of course not a proof. Much further work is 
required to support such a scenario with quantitative arguments — for example, it is not 
clear from this discussion that an attractive rather than a repulsive force will develop, even if 
Hlr does become unboundedly large for large separations between configurations of quarks 
that are not in color singlet configurations. Whether or not these scenarios for quark color 
confinement are borne out by further investigations, H — the QCD Hamiltonian in the M 
representation of the temporal gauge — promises to be a productive formulation of QCD 
dynamics for further investigations into the low-energy properties of QCD. 

This research was supported by the Department of Energy under Grant No. DE-FG02- 
92ER40716.00. 
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